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摘要 : 利用 分 析 法 研究 了 具有 磁 辐 射 的 两 成 分 模型 (党 层 和 中 子 超 流 体 ) 的 脉冲 星 在 可 
变 的 磁 辐 射 制 动 力矩 的 作用 下 ， 两 成 分 自 旋 角速度 随时 间 的 长 期 变化 。 给 出 了 具有 磁 辐 射 两 
成 分 模型 的 耦合 方程 组 的 分 析 解 。 理 论 结果 给 出 两 成 分 模型 在 外 力 可 变 的 磁 辐 射 制 动 力 抵 的 
作用 下 ， 自 旋 角 速度 随时 间 长 期 减 慢 。 利 用 所 得 的 分 析 解 对 具有 磁 辐 射 的 两 成 分 模型 蟹 状 星 
云 脉 冲 星 (PSR0531+21) ( Crab) 在 磁 辐 射 力 矩 可 变 的 情况 下 做 了 数值 计算 。 并 讨论 了 所 得 的 
理论 和 数值 结果 。 结 果 表 明 ， 艘 状 星云 脉冲 星 (PSR0531+21) 在 磁 辐 射 制 动 力矩 的 作用 下 ， 
索 层 自 旋 角速度 随时 间 的 长 期 减速 每 年 为 -0. 245 s。 
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为 了 解释 脉冲 星 突然 加 速 的 现象 ， 许 多 学 者 对 脉冲 星 提 出 各 式 各 样 的 模型 ， 其 中 较 成 功 的 模型 如 
磁 偶 极 辐射 模型 、 四 极 弹 性 能 模型 以 及 两 成 分 模型 等 ， 两 成 分 模型 最 先 由 文 [1-3] 提 出 。 两 成 分 模型 
主要 是 由 导电 较 高 的 固体 外 壳 和 壳 内 的 中 子 超 流体 的 混合 物 组 成 。 因 此 中 子 星 的 构造 模型 显然 分 为 两 
部 分 : 一 部 分 是 带电 的 固体 外 壳 ; 另 一 部 分 是 这 内 的 中 子 超 流体 的 混合 体 。 此 外 ， 两 成 分 以 不 同 角 速 
度 自 转 。 自 转 使 超 流 体 成 分 同 带 电 成 分 在 粘 汪 作用 下 相 耦 合 ， 尽 管 这 种 耦合 是 弱 的 ， 但 也 反应 出 超 流 
体 的 丰富 度 和 耦合 的 程度 。 由 于 脉冲 星 磁 偶 极 辐射 在 脉冲 星 表 层 产 生 磁 辐 射 制 动 力矩 ， 这 种 力矩 使 两 
成 分 自 旋 角速度 改变 。 因 此 ， 只 要 知道 磁 转 矩 变化 以 及 可 观测 的 宏观 驰 耶 时 间 + 和 中 子 超 流体 的 丰富 
度 QO， 两 成 分 的 自 旋 耦 合 的 微分 方程 可 解 。 文 [4] 给 出 了 对 两 成 分 耘 合 方程 组 在 磁 辐 射 制 动 力矩 不 变 
的 情况 下 耦合 方程 组 的 解 并 对 解 做 了 讨论 。 文 [5] 研 究 了 中 子 星 的 发 射 噪声 的 两 成 分 模型 , 但 没有 研 
究 磁 辐射 力矩 对 两 成 分 模型 角速度 改变 的 影响 。 文 [6] 人 研究 中 子 星 两 成 分 模型 在 广义 相对 论 构 染 内 的 
自转 动力 方程 ， 并 假定 过 成 分 的 转动 20. 为 常数 ， 而 中 子 超 流体 层 的 角速度 依赖 于 坐标 ， 但 没有 研究 
人 磁 辐 射 力矩 对 两 成 分 模型 角速度 改变 的 影响 。 本 文 作者 研究 了 在 磁 辐 射 制 动 力矩 随时 间 可 变 的 情况 
下 ,具有 人 磁 辐 射 的 两 成 分 角速度 变化 的 规律 ， 并 给 出 方程 组 的 解 。 最 后 将 解 应 用 于 具有 人 磁 辐 射 两 成 
分 模型 的 多 状 星云 脉冲 星 ( Crab pulsar) 的 两 成 分 角速度 随时 间 变 化 的 数值 解 ， 此 解 不 同 于 文 [4] 给 出 
的 解 。 


1 脉冲 星 两 成 分 的 耦合 方程 式 及 其 在 磁 辐 射 力矩 不 变 条 件 下 的 解 


根据 两 成 分 模型 理论 ， 两 成 分 的 自 旋 角速度 并 不 一 致 。 假 定 带电 的 固体 外 壳 的 转动 惯量 为 1 并 以 
角速度 Q2 绕 自转 轴 旋 转 ， 而 党 内 中 子 超 流体 的 转动 惯量 并 以 角速度 Q2, 绕 自转 轴 旋 转 , 但 NANA, 
带电 外 壳 成 分 与 过 内 中 子 超 流体 成 分 之 间 的 耦合 由 下 式 给 出 ” : 
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d0 I, 
B T aR (1) 
4, Loop 
hp cs ed . (2) 
(1) 式 、(2) 式 是 在 星 震 后 不 存在 跃 变 时 的 方程 组 。 其 中 ，N(b) 为 外 部 辐射 制 动 力矩 ; 7 是 理论 上 计 
算 的 脉冲 星 弛 豫 时 间或 称 微观 弛 驳 时 间 , 它 依赖 于 2。 文 [7] 给 出 计算 r. 的 公式 ， 文 [8] 将 r. 用 下 式 
Am: 
1 Qf/ A KT TA? 
T, 40 f E E, P (- ur) p (3) 


其 中 ，A 为 超 流体 间隙 参数 (能 量 间隔 ) ; Ej 为 电子 的 费 米 能 量 。 
微观 弛 豫 时 间 m. n T HZ ER IER] r 表示 ， 宏 观 弛 豫 时 间 > 是 可 以 观测 的 时 间 ， 两 者 之 间 关 系 


H FREER”, ; 
en »|l1 
u RII (4) 
> 中 子 超 流体 的 丰富 度 0: 
L AQ, 
e- i-a) (5) 
m 其 中 ，AQ, 为 0, 的 初始 跃 变 ， 而 总 转动 惯量 I: 
- I-I, +1.. (6) 
T 文 [8] 指 出 只 要 D,1,, Q 可 以 表示 : 
ps QsI/I. (7) 
© 由 于 大 多 数 脉 串 星 均 满足 上 述 条 件 ， 故 : 


N I, = QI, I 2I-1, 2I(1- Q). (8) 
> 因此 ，(4) 式 可 以 写成 ; 
二 = or， (9) 


X[4]28:H f TEREARET 012] N 不 变 的 情况 下 ， 耦 合 方程 组 (1) - (2) 的 解 : 


N L 
N=-7t t 0677 + (2,， 
(10) 
N 
0, -0- Q6 + 


以 下 本 文 给 出 在 辐射 制 动 力矩 可 变 的 情况 下 ， 方 程 组 (1) ~ (2) 的 解 并 应 用 于 脉冲 星 。 
2 可 变 磁 辐射 制 动 转 矩 的 形式 
本 文 研究 具有 磁 辐 射 的 两 成 分 模型 的 脉 串 星 。 有 磁 辐 射 必然 有 磁 辐 射 力矩 作 用 在 具有 了 磁 辐 射 的 两 


成 分 模型 上 ， 所 以 方程 组 (1) ~ (2) 中 N 是 磁 辐 射 力矩 。 本 文 给 出 磁 辐 射 力矩 随时 间 变 化 的 公式 。 
按 磁 偶 极 辐射 模型 ， 脉 冲 星 的 辐射 功率 W, 是 由 自转 能 的 变 率 E 转 化 来 的 ， 即 


W, +E =0 3 W, =-E. ao 
辐射 功率 和 自转 能 的 变 率 可 写成 ' "|. 
32m! WA ,dn 1 
z = E- IQ E = IČ. 12 
d 3c3P4 3c? , dt 9 2 ( ) 
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由 (11) 和 (12) 式 可 得 
dQ 2È, 
0 13 
"iE PI” i 


Ab, RIE u=R B sina; 尺 为 脉冲 星 半径 ; B, 为 脉冲 星 的 表面 磁场 。 假 定 磁 转 矩 垂直 于 自转 轴 ， 则 
ac=90 ， 故 磁 矩 w= 玉 B.。 如 果 在 短 时 间 内 不 考虑 磁 衰 减 或 增长 ， 则 (13 ) 式 成 为 
dO)  2P6( ， 


" 3 pe" 角 动 量 J 210. 
y dJ 2 7 3 
in une 14 
所 ~ T (14) 


JCGQ) T g 3 to (1c)? 

积分 式 的 下 限 如 EE ERE ANT EE EK E SE LAE R08 9] B PALLA SEA USB E TRE, AAWA 
IR] to 不 是 脉冲 星 诞生 开始 :=0 时 的 初始 时 间 ， 它 是 具有 磁 辐 射 的 两 成 分 模型 统一 体 从 现在 年 龄 开始 
的 时 间 。7(5) 是 二 时 的 角 动 量 : 


4 2p -1/2 
J)a =J 3). | , 
C 


i 2p p 
积分 上 式 : mM : 


3 P 
所 以 : JG), JODIE KG 78)]77, (15a) 
4Uj Aw 
其 - - 
其 中 ， Elus. t (15b) 
TI C15a) fICIT5b) nT A28 Hh RET E] 71x. NG) 的 表达 式 : 
N(t) =- Y, 
将 (15a) 式 代入 上 式 后 得 
N(t) Sed EE +K (t T t) ] ^ = 了 Ji [1 * K,(t E iud] 0. (16) 


3 — RET i 2] 7] XR n] 2E TED P PALLA 21 REEH CL) ~ (2) 的 分 析 解 


将 (8) (9) RI(16) 5X9 1, ,. 1, 和 E 以 及 NG)IAERRÉBO) ~ (2) 式 ， 则 (1) ~ (2) 式 可 写成 如 
FÉR: : 


dO 1 JK E de. Qp 
RI (17) 
dQ Š 
a A y EA (18) 
dt T 
这 样 ， 方 程 组 (1) ~ (2) nT AIE EBRE] 7 表示 。 将 方程 组 (1) 和 (2) 的 两 端 各 项 相 加 ， 可 得 
dQ , dO, d M 
Lowe. a 04 + LA) -- N(t). 


转动 惯量 I 和 工 不 变 或 为 常数 的 情况 下 ， 将 (16) 式 的 N(t) 代 入 上 式 积分 后 可 得 


t 1 t 
40,0 + 1,0) 5- [voy =- SKJ (00) [1 KG ~ 177, 


to 


所 以 : JO), ELOI, -J(4),L1-K(t-1)] ^. (19) 


280 天 文 研究 与 技术 13 卷 


此 式 正好 是 (15a) 式 ， 这 说 明 角 动量 变化 的 (19) 式 左 端的 (1), 是 两 成 分 角 动 量 之 和 。 
利用 (7) ~(8) 式 将 (19) 改 成 下 列 两 式 : 


_J(to) aan h E JC) RS -M2 . Q . 
N= I [1 +K,(t-t)]7 ra- t5 i2] (29 9s (20) 
je moglesqii eae abren a (21) 


1, Q Q 
3b, 10) 20), | BG). 10 (0. 


将 (21) 式 的 Q 代入 (17) 式 ， 方 程 组 (17) ~ (18) 变 成 下 列 一 组 一 阶 非 线性 齐 次 方程 组 
dQ 


dt 


+ AA) [1 -tK(t-4)]"?-OQ(nu),— 


3 agU *tK(t-4)] ^, (22) 


ITIN C 
dt T i T 
解 (22) 式 得 到 Q 后 再 代入 (23) 式 得 到 0 的 解 。 
同样 ,将 (20) 的 Q 代 入 (18) 式 , 方程 组 (17) ~ (18) 式 变 成 下 列 另 一 种 形式 : 


0. (23) 


dQ, 1 O (to) m 

t= = [1 +K (re$). (24) 
dt T T 

dO Q,, 1 Q(&),K, a2, 

d e 7g ugy (DNUS CRUUIPC A olds 25) 


方程 组 (22) ~ (23) 同方 程 组 (24) ~ (25) 的 解 是 等 价 的 。 以 下 选取 方程 组 (22) - (23) RE LQ 和 0Q,。 首 
先 解 一 阶 线性 非 齐 次 方程 (22) : 


+ PQ = Wa). 
POEL, WO) 206), 0 RG tk) 
7 20 - 9) 
解 的 形式 为 
a-«F*[acp.[| Lt o1 ^- Tg AU to) rta e). 
4. A=, +L) + Cew™. (26) 


积分 式 : 


1 
1, 2 Q0, SUL + KG o 8)] ed 


= Aio) e^ [i + KG 70377 * KE KG 78)17? * T GSD?U KG 78)]177] , 


L=- Q(t), 1 +K, (t-t) ] ^e" dt 


K, 
2(1- ii -jll 


hato t [1 +K- 19102 4.2. O80" 
2(1 - Q) ) 4 (1-09) 


I, f LIA (26) X, 
N= ALAK G-t) ^ -EK(L*K(r-4)] ^ -K;i1-K(r-5)] ^|] * Ce", (27) 


=- (X(t,) mel 


[1 + 天 (一 add 
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pr. TTE ES TA 

cil 2 1-9 4 1-9Q 
利用 初始 条 件 t=t0, Q1) 2 2 (to),,， 这 表示 磁 辐 射 和 两 成 分 统一 体 在 星 震 和 跃 变 后 从 现在 辐射 年 龄 
to 开始 , 则 由 上 式 可 得 


4 
Hp, 人 = 可 


(28) 


C, =A) (K, + Ke. 
所 以 壳 层 角速度 2 在 磁 辐 射 力矩 作用 下 随时 间 变 化 是 
Q-20()([1*K(r-1)] ^ -Kjy1*K(t-1)] ^? 
E K,[1 +K, (t i ae * (K, *Kje C7") . (29) 


EN = 0(1) - (1). (30) 
TE (29) XP, i01; m (6 8AD) H7 At, At 是 以 年 为 单位 的 时 间 间 隔 , 它 与 脉冲 星 年 龄 无 关 , 也 与 起 始 时 
间 无 关 ， 所 以 (29) 式 可 以 写成 用 时 间 间 隔 表 示 的 式 子 : 
Q-2O0(4)([1* K(A] ^ - K/y|1T- K,At] °? -K[1 +K At] + (K, t Ky)e 77) , (31) 
此 式 只 与 取 的 时 间 间 隔 有 关 ， 而 与 时 间 起 点 无 关 。 故 本 文 在 计算 时 取 时 间 间 隔 为 一 年 。 
将 (29) 式 代入 (23) 式 后 可 得 到 O, 的 一 阶 线性 非 齐 次 方程 , 令 其 解 为 
0, mg MEDIE 十 L, + L + L] FS Qu Umm (32) 
TI FIR TRAE ET D. Lu Lu LARMER EA C, ERUAVOA EX, HR ETUR 
流体 的 角 速 2 随时 间 变 化 的 解 (推导 式 子 较 长 ， 因 篇 幅 关 系 略 去 ) : 


0,2 0, [L1 * Ka 70177 A p Kerl + KG 710177 (R K) er 
e in ey Kr T ERO m qeu £ Q (1), e 070077 (33) 
60, =N), - Q(t), . (34) 


(29) 式 和 (33 ) 式 就 是 本 文 给 出 的 脉冲 星 磁 辐射 力矩 对 两 成 分 模型 自转 角 速 随时 间 长 期 变化 的 式 子 。 


4 脉冲 星 PSR0531+21( Crab) 的 理论 数值 结 


本 文 研究 的 脉冲 星 必须 具有 磁 辐 射 的 两 成 分 模型 ， 根 据 文 [ 3] ， 和 你 状 星云 脉冲 星 ( Crab ，PSR0531 
+21) 是 两 成 分 模型 ， 根据 文 [4] 它 又 是 磁 辐 射 模型 或 是 具有 这 两 种 模型 的 统一 体 。 本 文选 取 Crab 
脉冲 星 (PSR0531+21) 作 为 计算 实例 。 首 先 给 出 这 个 脉冲 星 的 物理 参数 如 表 1， 其 中 QO,、Q、7 引 自 
X[3, 10], ， 表 面 磁 场 B. 引 自 文 [11], j=R'B.。 一 般 假定 中 子 星 的 半径 R=1.2 x 10' cm， 转 动 惯量 
121.4 x 10?( em? * g) 


表 1 和 蟹 状 星云 脉冲 星 ( PSR0531+21) 的 物理 数据 
Table 1 Data for Pulsars PSR0S31+21 


脉冲 星 Pulsar N (rad/s) Q-I,/I T(d, yr) B(G) MX1030(cm * G) 


PSROS31+21( Crab) 190 0. 96 7.7 (d) 0. 999 553 952E+12 6. 53 


将 表 1 中 的 数据 代入 (28) 式 ， 得 到 PSR0531+21 BB K, =5.43 x 10 , K, 24.33 x 10^, K, =2.35 x 10°, 
再 取 时 间 间 隔 Ar=1 年 。 将 有 Ka, K ZURIA (30) (31) 3,442] PSR0531+21 (Crab) 的 壳 层 在 
磁 辐 射 力矩 作用 下 角 速 每 年 变化 的 数值 如 表 2。 
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R2 和 蟹 状 星 云 脉 冲 星 (PSR0531+21) 的 数值 的 结果 (每 年 变化 的 数值 : At=1 5E) 
Table 2 The numerical results for PSR0531+21 per year ( At=1 year) 


脉冲 星 Pulsar (1t) (rad/s) AYD 602( rad/s) 


PSR0531+21( Crab) 189. 755 0 0. 999 871 -0. 245 0 
人 2 的 数值 在 表 1 中 给 出 。 


对 于 中 子 超 流体 的 角速度 变化 情形 由 于 初始 角速度 Cu) ,难以 观测 到 , 所 以 中 子 超 流体 的 角速度 
ACE) ,随时 间 变 化 也 同样 难以 观测 到 。 故 本 文 只 能 根据 所 推 得 的 理论 式 子 对 壳 层 角速度 变化 做 一 计 
算 ， 但 给 出 中 子 超 流体 的 角速度 变化 的 理论 式 子 仍 有 理论 价值 和 意义 。 


5 讨论 和 结论 


(1) 对 推 得 的 理论 结果 式 子 的 验证 。 将 初始 时 间 t=t 代 人 理论 式 子 (29) 和 (33) ， 可 得 到 右 端的 
初始 角速度 Q(t) 和 Q(t6),。 即 (29) ~(30) 和 (33) ~ (342) RER 2 0(1,) ,602=0;0 =A (t), ôN, = 
0。 这 说 明文 中 给 理论 结果 是 正确 的 。 此 外 ，(29) 和 (33) 式 的 级 数 随时 间 增 大 是 收敛 的 。 

(2) 根 据 表 2 给 出 的 数值 结果 ， 脉 冲 星 PSR0531+21 在 磁 辐 射 外 力矩 的 作用 下 ， 外 壳 角 速度 随时 
间 逐 渐 减 慢 ， 然 而 这 种 减速 是 长 期 的 ， 而 和 由 于 壳 震 角 速 突 然 加 速 的 路 变 两 者 过 然 不 同 。 前 者 是 长 期 
性 的 ， 后 者 是 突然 临时 性 的 。 此 外 ， 前 者 长 期 减速 不 会 影响 后 者 突然 跃 变 的 加 速 ， 而 后 者 的 突然 跃 变 
加 速 对 前 者 的 长 期 变化 也 不 产生 影响 。 因 为 方程 组 (1) ~ (2) 是 在 星 震 跃 变 不 存在 的 情况 下 成 立 ， 
此 它 的 解 (长 期 变化 ) 不 受 星 震 和 跃 变 的 影响 。 

(3) 脉冲 星 磁 辐 射 的 演化 起 点 是 从 脉冲 星 诞生 := 0 开始 ， 而 两 成 分 模型 的 演化 起 点 是 从 星 震 和 跃 
变 后 为 起 点 ,两 者 并 不 一 致 。 本 文 致力 于 如 何 使 两 个 模型 的 演化 起 点 合 二 为 一 ， 这 是 本 文 解决 此 问题 
的 特点 。 本 文 研究 的 脉冲 星 具 有 磁 辐 射 和 两 成 分 的 统一 体 ( Crab. 脉冲 星 就 是 这 两 种 模型 的 统一 体 ) , 
因此 ， 演 化 的 起 点 必须 在 星 震 和 路 变 后 统一 体 的 初始 时 间 。 如 果 按 磁 辐 射 的 初始 时 间 是 脉冲 星 诞 
^E t=0 为 起 始 时 间 ， 可 是 脉冲 星 诞 生 时 的 物理 参量 是 一 个 不 确定 的 较 大 的 物理 量 。 另 外 也 不 符合 两 
成 分 模型 的 演化 起 点 。 因 此 应 该 采用 脉冲 星 的 现在 年 龄 i, 为 两 个 模型 统一 体 的 初始 时 间 。 在 积分 式 子 
(29)~(30) 中 ， 取 从 六 到 1 为 积分 上 下 限 , 可 是 1=t6+At， 故 在 (29) ~ (30) 式 子 中 取 :0 2 (AC) -t= 
At。 按 本 文 计算 取 时 间 间 隔 At=1 年， 所 以 演化 时 间 间 隔 与 脉冲 星 演化 的 起 始 时 间 无 关 。 这 样 ， 具 有 
磁 辐 射 的 两 成 分 模型 ， 如 Crab 脉冲 星 的 演化 时 间 可 选取 在 星 震 和 跃 变 后 磁 辐 射 和 两 成 分 的 统一 体 的 
现在 年 龄 为 两 者 的 初始 时 间 。 这 个 方法 解决 了 上 述 两 个 模型 演化 起 点 不 一 致 的 问题 。 

(4) 本 文 假 定 (13) 式 中 的 磁 偶 极 矩 为 常量 ， 即 不 随时 间 误 减 或 增长 。 实 际 上 ， 如 果 从 长 期 考虑 
REB u 是 随时 间 变 化 的 。 即 凡 =Hoe,E=2《rn。 文 [9] 给 出 的 rr=1.6x10 年 ,上 =1.3x10 年 。 因 此 ， 
磁 误 减 时 间 较 长 ， 磁 衰减 系数 每 年 很 小 。 本 文 只 计算 每 年 角 速 的 变化 值 ， 因 此 ， 磁 衰减 的 影响 可 以 不 

男 一 方面 ， 有 的 脉冲 星 的 磁场 是 增长 的 。 例 如 文 [12] 给 出 的 PSR0531+21( Crab 脉冲 星 ) JE FR 
场 增长 的 脉冲 星 ， 即 久 =joe** =RiBoe”*， 但 磁场 增 量 很 小 ， 增 量 为 +0. 000 6 x 1012 G/ 年 。 这 样 小 的 
增 量 对 Crab 脉冲 星 每 年 壳 层 角速度 变化 的 影响 也 可 不 必 考 虑 。 

(5) 在 第 3 节 中 厂 、 厂 、 厂 大 各 积分 式 子 利用 了 下 面 的 积分 公式 : 

[ied n xii BUR D e — i 


a a 


级 数 展开 式 应 用 于 本 文 n=12, a=1/Kir Mx=[1+K(t-t)]o X 1P K=5.43 x 10 , r=7.7 d=665 280 s, 
K,T=3.6124Xx10°, x>1 随时 间 4 延长 而 增加 。 将 这 些 代 入 上 式 , 此 级 数 是 收 剑 的。 所以， h 
1;、 各 积分 式 是 收敛 的 ， 故 此 所 得 (29) 式 也 是 收敛 的 。 男 外 ， 如 果 将 第 4 节 表 1 FARK, K, K 
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的 数值 代入 (29) 式 也 可 得 到 (29) 式 是 收敛 的 ， 当 [1+Ki(t-io)] 随 时 间 ; 延 长 而 增 大 。 但 由 于 在 级 数 
中 上 略 去 了 [1+K,(t-t6)] “以 上 的 高 阶 项 ， 所 以 , 解 (29) 式 是 近似 的 。 
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The Impact of Magnetic Radiation Braking Torque on the Secular 


Retardation of Spin of Two-Components of Pulsar 


Li Linsen 
(School of Physics, Northeast Normal University, Changchun 130024, China, Email; dbsd-lls@ 163.com) 


Abstract: The secular influences of the variable magnetic radiation braking torque on the spin down of 


two-components of pulsar are studied by using analytical method. The analytical solutions of the coupled 


equation system are given with due consideration of the time-dependent magnetic braking torque by solving the 


system of non-homogeneous equations. The analytical solutions obtained are applied to the research for 


PSROS31+21 (Crab Pulsar). The numerical results are obtained in Table 2. The numerical results show that 


the spin of Crab pulsar speeds down -0. 245rad/yr. Discussions are held on the obtained results. 
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